Quantum fluctuations of systems of interacting electrons in two spatial dimensions. 
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1 Abstract 

The random matrix ensembles (RME) of quantum statistical Hamiltonian operators, e.g. 
Gaussian random matrix ensembles (GRME) and Ginibre random matrix ensembles (Gini- 
bre RME), are applied to following quantum statistical systems: nuclear systems, molecular 
systems, and two-dimensional electron systems (Wigner-Dyson electrostatic analogy). Mea- 
sures of quantum chaos and quantum integrability with respect to eigenergies of quantum 
systems are defined and calculated. Quantum statistical information functional is defined 
as negentropy (opposite of entropy or minus entropy). The distribution function for the 
random matrix ensembles is derived from the maximum entropy principle. 

2 Introduction 

Random Matrix Theory (RMT) studies as an example random matrix variables H, corre- 
sponding to random quantum Hamiltonian operators IK. Their matrix elements 7iij,i,j = 
1,...,N,N > 1, are independent random scalar variables and are associated to matrix ele- 
ments Hij of Hamiltonian operator H. [1, 2, 3, 4, 5, 6, 7, 8]. We will use the following nota- 
tion: H is hermitean Hamiltonian operator, H is its hermitean matrix representation, E^s 
are the eigenvalues of both H and H. The matrix elements of if are Hij, i, j — 1, N, N > 1. 
Moreover, K is random hermitean Hamiltonian operator, K is its hermitean random matrix 
representation. The values assumed by K are denoted by H, whereas the values of K are 
H , respectively. The SiS are the random eigenvalues of both K and K. The random matrix 
elements of K are Tiij, i,j = 1, N, N > 1. Both Si and TCij are random variables, whereas 
Ti, is random matrix variable. The values assumed by random variables Si are denoted Ei. 
There were studied among others the following Gaussian Random Matrix ensembles GRME: 
orthogonal GOE, unitary GUE, symplectic GSE, as well as circular ensembles: orthogo- 
nal COE, unitary CUE, and symplectic CSE. The choice of ensemble is based on quantum 
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symmetries ascribed to the Hamiltonian operator H. The Hamiltonian operator H acts on 
quantum space V of eigenfunctions. It is assumed that V is iV-dimensional Hilbert space 
V = F N , where the real, complex, or quaternion field F = R, C,H, corresponds to GOE, 
GUE, or GSE, respectively. If the Hamiltonian matrix H is hermitean H = W, then the 
probability density function fa of TC reads: 

f H (H)=C n(3 ex V [-P-yTT(H% (1) 

C W = {^f^\ 

M W = N+^N{N-1)P, 



J f n (H)dH 



N N D-l 

i=l j>i 7=0 

H ij = (H§\...,H^- 1) )eF 1 

where the parameter (3 assume values (3 = 1,2,4, for GOE(iV), GUE(iV), GSE(A r ), respec- 
tively, and J\fn/3 is number of independent matrix elements Hij of hermitean Hamiltonian 
matrix H. The Hamiltonian H belongs to Lie group of hermitean N x N dimensional F- 
matrices, and the matrix Haar's measure dH is invariant under transformations from the 
unitary group U(iV, F). The eigenenergies Ei,i = 1,...,JV, of H, are real-valued random 
variables Ei = E*, and also for the random eigenenergies it holds: £j = £*. It was Eugene 
Wigner who firstly dealt with eigenenergy level repulsion phenomenon studying nuclear 
spectra [1, 2, 3]. RMT is applicable now in many branches of physics: nuclear physics (slow 
neutron resonances, highly excited complex nuclei), condensed phase physics (fine metallic 
particles, random Ising model [spin glasses]), quantum chaos (quantum billiards, quantum 
dots), disordered mesoscopic systems (transport phenomena), quantum chromodynamics, 
two-dimensional Euclidean quantum gravity (2D EQG), Euclidean field theory (EFT). 



3 The Ginibre ensembles 

Jean Ginibre considered another example of GRME dropping the assumption of hermiticity 
of Hamiltonians thus defining generic F-valued Hamiltonian matrix K [1, 2, 9, 10]. Hence, K 
belong to general linear Lie group GL(iV, F), and the matrix Haar's measure dK is invariant 
under transformations from that group. The distribution fa of /C is given by: 

f K (K) = C^exp [-P-y Tt(K^K)], (2) 
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Mkp = N 2 f3, 
J f K {K)dK = 1, 

N N D-l 
i=l j=l 7=0 

where (3 = 1,2,4, stands for real, complex, and quaternion Ginibre ensembles, respectively, 
and JC is random matrix variable associated with matrix K. Therefore, the eigenenergies Zi 
of quantum system ascribed to Ginibre ensemble are complex-valued random variables. The 
eigenenergies Zi,i = 1,...,N, of nonhermitean random matrix Hamiltonian /C are not real- 
valued random variables Z { ^ Z*. Jean Ginibre postulated the following joint probability 
density function of random vector Z of complex eigenvalues Zi, Z^ for Nx N dimensional 
random matrix Hamiltonians /C for (3 = 2 [1,2, 9, 10]: 

P(Z 1 ,...,Z N )= (3) 

N -i N N 

= n^-n^-^f •ex P (-£i^i 2 ), 
j=i j- i<j j=i 

where Zi are complex- valued sample points (values) of Z; L {Z i e C). 

We emphasize here Wigner and Dyson's electrostatic analogy. A Coulomb gas of iV unit 
charges Qi moving on complex plane (Gauss's plane) C is considered. The complex vectors 
of positions of charges are Zi and potential energy U of the system is: 

N i 

U(Z ± , Z N ) = - £ In \Zi - Z 3 \ + - £ \Z t \ 2 . (4) 

i<j i 

If gas is in thermodynamical equilibrium at temperature T = {(3 = = 2, k B is Boltz- 
mann's constant), then probability density function of vectors Zi of positions is P(Zi, Z^) 
Eq. (3). Therefore, complex eigenenergies Zi of quantum system are analogous to vectors of 
positions of charges of Coulomb gas. Moreover, complex- valued spacings A 1 ^ (first order 
forward/progressive finite differences) of complex eigenenergies Zi of quantum system: 

A 1 ^ = AZi — Z i+1 — Zi,i — 1, (N - 1), (5) 

are analogous to vectors of relative positions of electric charges. Finally, complex- valued sec- 
ond differences A 2 Zi (second order forward/progressive finite differences) of complex eigenen- 
ergies Zf 

A 2 Z, t = Z t+2 - 2Z i+1 + Z, t ,i = 1, (N - 2), (6) 

are analogous to vectors of relative positions of vectors of relative positions of electric charges. 

The eigenenergies Zi = Z{%) can be treated as values of function Z of discrete parameter 
% — 1, ...,N. The "Jacobian" of Zi reads: 

JacZ (i , = Ja cZ, = f = § = ^§ = A% (7, 
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where Ai — i + 1 — 1 = 1. We readily have, that the spacing is an discrete analog of Jacobian, 
since the indexing parameter i belongs to discrete space of indices i & I — {1,...,N}. 
Therefore, the first derivative ^ i with respect to % reduces to the first forward (progressive) 
differential quotient The Hessian is a Jacobian applied to Jacobian. We immediately 
have the formula for discrete "Hessian" for the eigenenergies 

He»Z(i) = He SS Z, = = K = g| = AX (8) 

Thus, the second difference of Z is discrete analog of Hessian of Z. One emphasizes that both 
"Jacobian" and "Hessian" work on discrete index space / of indices i. The spacing is also a 
discrete analog of energy slope whereas the second difference corresponds to energy curvature 
with respect to external parameter A describing parametric "evolution" of energy levels 
[11, 12]. The finite differences of order higher than two are discrete analogs of compositions 
of "Jacobians" with "Hessians" of Z. 

The eigenenergies E i} i e /, of the hermitean Hamiltonian matrix H are ordered in- 
creasingly real numbers. They are values of discrete function = E{i). The first order 
progressive differences of adjacent eigenenergies: 

A 1 E i — E i+1 — Ei,i — 1, (N - 1), (9) 

are analogous to vectors of relative positions of electric charges of one-dimensional Coulomb 
gas. It is simply the spacing of two adjacent energies. Real- valued progressive finite second 
differences A 2 Ei of eigenenergies: 

A 2 Ei = E i+2 - 2E i+l +E h i = 1, (N - 2), (10) 

are analogous to vectors of relative positions of vectors of relative positions of charges of 
one-dimensional Coulomb gas. The A 2 Zi have their real parts ReA 2 Zj, and imaginary parts 
ImA 2 Zj, as well as radii (moduli) \A 2 Zi\, and main arguments (angles) ArgA 2 Z«. A 2 Zi are 
extensions of real-valued second differences: 

Re(A 2 Z,) = Re(Z l+2 - 2Z, l+1 + Z % ) = A 2 ReZ t = A 2 E t ,i = 1, .., (N - 2), (11) 

of adjacent ordered increasingly real- valued eigenenergies of Hamiltonian matrix H defined 
for GOE, GUE, GSE, and Poisson ensemble PE (where Poisson ensemble is composed of 
uncorrelated randomly distributed eigenenergies) [13, 14, 15, 16, 17, 18, 19, 20, 21, 22]. The 
"Jacobian" and "Hessian" operators of energy function E(i) = Ei for these ensembles read: 

t jpf\ t T7i dEi dE % A 1 ^ AEi . 

JacE W = J acEl = — = — * — = — = AEk, (12) 



and 



HessE(,) = Ress Ei = ^ = ^- ^ = W = A ^- ( 13 ) 
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The treatment of first and second differences of eigenenergies as discrete analogs of Jacobians 
and Hessians allows one to consider these eigenenergies as a magnitudes with statistical 
properties studied in discrete space of indices. The labelling index i of the eigenenergies 
is an additional variable of "motion", hence the space of indices / augments the space of 
dynamics of random magnitudes. 

One may also study the finite expressions of random eigenenergies and their distributions. 
The finite expressions are more general than finite difference quotients and they represent 
the derivatives of eigenenergies with respect to labelling index i more accurately [23, 24]. 

4 The Maximum Entropy Principle 

In order to derive the probability distribution in matrix space M. = MATRIX (iV, N, F) of 
all N x N dimensional F-valued matrices X we apply the maximum entropy principle: 

max{S p \f x ]:{l) = l,(h x ) = Up}, (14) 

whereas (...) denotes the random matrix ensemble average, 

(g x ) = I g x (X)f x (X)dX, (15) 

and 

Sp[f x ] = (-k B lnf x ) = [ (-k B lnf x (X))f x (X)dX, (16) 
is the entropy functional, 

(h x ) = I h x (X)f x (X)dX, (17) 

and 

(1) = / lf x (X)dX = I. (18) 

Here, h x stands for "microscopic potential energy" of random matrix variable X, and Up is 
"macroscopic potential energy". We recover the Gaussian random matrix ensemble distri- 
bution fa Eq. (1), or Ginibre random matrix ensemble distribution fa Eq. (2), respectively, 
if we put h x (X) = Tr(^X^X), and X = H, or X = K, respectively. The maximization of 
entropy Sp[f x ] under two additional constraints of normalization of the probability density 
function f x of X, and of equality of its first momentum and "macroscopic potential energy", 
is equivalent to the minimization of the following functional T\j x \ with the use of Lagrange 
multipliers oi\,{i\. 

min{^[/^]}, (19) 
F[fx]= [ (+k B \nf x (X))f x (X)dX + a 1 [ f x (X)dX 

J M J M 

+A J h x (x)f x (x)dx = 

= -Sp[fx] + ai (1> + A (h x ). 
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It follows, that the first variational derivative of T\fx\ must vanish: 

^1=0, (20) 

which produces: 

k B (hif x (X) + l) + a 1 + l3 1 h x (X) = 0, (21) 

and equivalently: 

f x [X]=C X f,-exp[-l3-h x (X)] (22) 
Cxa = exp[-(ai + 1) • k B 1 ],p = fa • k B \ 

One easily shows that 



6,1 >a (23) 
The variational principle of maximum entropy does not force additional condition on func- 
tional form of h x (X). The quantum statistical information functional Ip is the opposite of 
entropy: 

Wx] = Sp[fx\ = [ (+k B \nf x (X))f x (X)dX. (24) 

Information is negentropy, and entropy is neginformation. The maximum entropy principle 
is equivalent to minimum information principle. 
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